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Central  Limit  Theorems  in  the  Area  of  Large 
Deviations  for  Some  Dependent  Random  Variables 


by 

Narasinga  Rao  Chaganty  and  Jayaram  Sethuraman 


ABSTRACT 


A  magnetic  body  can  be  considered  to  consist  of  n  sites,  where  n  is  large. 

The  magnetic  spins  at  these  n  sites,  whose  sum  is  the  total  magnetism  present  in 

the  body,  can  be  modelled  by  a  triangular  array  of  random  variables  (x|n^ , . . . ,X^n^) 

Standard  theory  of  Physics  would  dictate  that  the  joint  distribution  of  the  spins 

can  be  modelled  by  dQ^x)  =  z"1  exp[-Hn(x)]  IIdP(x..)  ,  where  x  =  (Xj,...,xn)  e  Rn, 

where  Hr  is  the  Hamiltonian,  zn  is  a  normalizing  constant,  and  P  is  a  probability 

measure  on  R.  For  certain  forms  of  the  Hamiltonian  H  ,  Ellis  and  Newman 

n 

(Z.  Wahrscheinlichkeitsheorie  und  Verw.  Gebiete  44(1978)  117-139)  showed  that 

1-1/2 

under  appropriate  conditions  on  P,  there  exists  an  integer  rsl  such  that  Sn/n 

converges  in  distribution  to  a  random  variable.  This  limiting  random  variable 

is  Gaussian  if  r=l  and  non-Gaussian  if  r£2.  In  this  article,  utilizing  the 

large  deviation  local  limit  theorems  for  arbitrary  sequences  of  random  variables 

of  Chaganty  and  Sethuraman  (Ann.  of  Probability,  13  (1985)),  we  obtain  similar 

central  limit  theorems  for  a  wider  class  of  Hamiltonians  H  ,  which  are  functions 

n 

of  moment  generating  functions  of  suitable  random  variables.  We  also  present 
a  number  of  examples  to  illustrate  our  theorems. 


1.  Introduction. 

-  In  this  article  we  obtain  central  limit  theorems  for  some  dependent  random 
variables  which  are  used  to  describe  the  distribution  of  magnetic  spins  present 
in  a  ferromagnet  crystal.  A  ferromagnet  crystal  consists  of  a  large  number  of 
sites.  The  amount  of  magnetic  spin  present  at  site  i  will  be  denoted  by 
X{n) ,  i  =  1,  ...»  n,  where  n  is  a  positive  integer.  The  magnetic  spin  present 
at  any  site  interacts  with  the  magnetic  spins  at  its  neighboring  sites  and 
hence  gives  rise  to  some  dependency  among  the  random  variables  X^'s.  In  the 
Ising  model,  the  joint  distribution,  at  a  fixed  temperature  T  >  0,  of  the  spin 
random  variables  (X^  , . . .  ,X^n^ ) ,  is  given  by 

r  h  mi 

/  i  i  \  irt  /■  \  1  I  ^  I  tt  \ 


(1.1)  dQ^x)  =  exp  I- — — J  ndp(xj)> 

where  x  =  (x^,...,xn)  e  R°  and  P  is  a  probability  measure  on  R.  The  function 
Hn(x)  is  known  as  the  Hamiltonian  and  it  represents  the  energy  of  the  crystal 
at  the  configuration  x,  and  zn  is  a  normalizing  constant  which  is  also  known  as 
the  partition  function.  In  many  cases,  an  explicit  evaluation  of  zr  is  very 
difficult  and  physicists  usually  try  to  evaluate  the  limiting  free  spin  per  state 
£(T) ,  at  the  temperature  T,  defined  as  follows: 


(1.2) 


C(T)  =  -  li™  [(log(zn))/n3. 


For  some  particular  types  of  Hamiltonians,  it  has  been  shown  by  physicists  that 
there  exists  a  temperature  level  Tc  such  that  the  function  £(T)  is  infinite  or 
finite  according  as  T  <  T^  or  T^Tc  (see  Kac  (1968)).  A  phase  transition  is  said 
to  occur  at  the  critical  temperature  T^.  As  pointed  out  by  Ellis  and  Newman  (1978b), 
the  existence  of  the  critical  temperature  can  be  demonstrated  in  yet  another  way. 

We  may  be  able  to  show  that  for  T  >  T^,  there  is  a  weak  dependence  among  the  random 
variables  (X^  , . . .  ,X^n^ )  and  a  standard  central  limit  theorem  is  valid  for  Sn/>^n 
and  that  for  T  =  T^.,  there  exists  a  6  e  (1,2)  such  that  Sn/n^^  converges  to  a 


A  > 


't  l  ■ 
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non-Gaussian  limit  and  for  T  <  T^,  due  to  the  strong  dependence  among  the  xP^'s, 
the  random  variables  tend  to  cluster  in  several  ergodic  components.  This  shows  a 
marked  discontinuity  in  the  asymptotic  distribution  of  Sn  as  the  temperature  T  is 
allowed  to  vary  and  represents  our  approach  to  demonstrating  a  phase  transition. 

In  section  2,  we  consider  a  special  case  for  the  Hamiltonian,  Hn,  by  setting 
it  to  be  equal  to  -  EEx^Xj .  This  is  known  as  the  Curie-Weiss  model.  The 
asymptotic  distribution  of  Sn  for  this  model  when  P  (which  appears  in  Theorem  2.1) 
is  symmetric  Bernoulli  is  obtained  by  Simon  and  Griffiths  (1973) .  In  a  two  paper 
series,  Ellis  and  Newman  (1977,  1978a,  1978b)  extended  Theorem  2.1  of  Simon  and 
Griffiths  to  the  class  of  probability  measures  L,  defined  in  2.2  (see  also  Ellis 
and  Rosen  (1979)).  We  state  their  extension  precisely  in  Theorem  2.6.  Recently 
Jong-Woo  Jeon  (1979)  in  his  Ph.D.  dissertation  gave  a  simpler  and  statistically 
motivated  proof  of  Theorem  2.6.  The  goal  of  this  article  is  to  extend  Theorem  2.6 
for  a  larger  class  of  Hamiltonians  Hn  and  probability  measures  P.  Our  main  result. 
Theorem  3.7  is  stated  in  Section  3.  The  proof  of  Theorem  3.7  relies  on  a  recent 
large  deviation  local  limit  theorem  of  Chaganty  and  Sethuraman  (1985) ,  which  is 
restated  in  Section  3  as  Theorem  3.4  . 

Let  Tn,  n  2  l,  be  an  arbitrary  sequence  of  random  variables  with  analytic 
moment  generating  function  $n(z).  We  assume  that  Tn  satisfies  the  conditions  of 
Theorem  3.4  .  In  our  generalized  model  the  Hamiltonian  Hn(x)  is  taken  to  be 
equal  to  -  log(4>n(sn/n)],  where  s^  =  x^  +  . . .  +  xn>  Thus,  the  joint  distribution 
of  the  spin  random  variables  (X^  , . . .  ,X^n^)  is  given  by 

(1.3)  dQ^x)  =  2nX  <*>n(-sn/n)  ndpCxj)  » 

where  P  is  an  arbitrary  probability  measure.  Let  Sn  =  X^+  ...  +  X^n^ .  Under 
appropriate  conditions  on  the  probability  measure  P  we  show  in  Theorems  3.7  and 
3.18,  there  exists  an  integer  r  £  1  such  that  Sn/n*”^r  converges  in  distribution 
to  a  random  variable  Yr* ,  which  has  a  non-normal  distribution  when  r  s  2  and  normal 
distribution  when  r  =  l.  The  technique  of  our  proof  is  to  introduce  a  new  random 
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variable  W  ,  conditional  on  which,  X^n\  become  i.i.d.  It  is  easy 

to  obtain  the  limiting  distribution  of  WR  and  the  conditional  asymptotic  distri¬ 
bution  of  Sn/n*~*^r.  Using  the  results  of  Sethuraman  (1961)  we  deduce  the 
asymptotic  distribution  of  Sn/n^~^r. 

We  now  briefly  give  our  reasons  for  calling  these  theorems  on  the  asymptotic 
distribution  of  Sn  under  Q^,  defined  in  (1.3),  as  limit  theorems  in  the  area 
of  large  deviations.  A  standard  technique  to  obtain  the  asymptotic  distribution 
of  Sn  under  0  is  to  first  obtain  the  asymptotic  distribution  of  Sn  under  Pn, 
where 

(1.4)  dPn(x)  =  ndP(x.) 

and  then  to  use  contiguity  arguments,  as  in  LeCam  (1960).  This  technique  breaks 

down  completely  if  rs2.  For  the  various  models  considered  in  Physics  which 

are  described  in  greater  detail  in  Sections  2  and  3, 

d0  (x)  H  (x) 

Cl. 5)  |Ln(x)|  =  I  log  d?TxT  1  =  1  -V11-  -  lo*  znl 

converges  to  “  in  probability  under  Pfl  and  thus  contiguity  arguments  are  not 
applicable  here.  Under  Pn,  Sn/i/n  has  a  limiting  normal  distribution,  also, 
|Ln(x)|  is  small  in  the  area  of  ordinary  deviations  of  Sn,  that  is,  when  Sn/VrT 
is  finite,  while  it  is  large  otherwise.  Thus  from  the  point  of  view  of  P  ,  we 
are  looking  for  the  asymptotic  distribution  of  Sn,  when  Pn  is  modified  by  Lr(x) , 
which  is  substantially  different  from  1  in  the  area  of  large  deviations  of  S^. 
This  view  point  helps  in  a  statistically  motivated  proof  of  the  asymptotic 
distribution  of  under  and  describes  the  background  behind  the  title  of 
this  article.  One  should  also  note  that  the  normalizing  factor  on  in  its 
asymptotic  distribution  under  0  is  different  from  the  corresponding  factor 


under  P  . 
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2 .  A  Brief  Summary  of  Curie-Weiss  Model  and  Its  Extensions. 

In  a  ferromagnetic  system  with  only  isotropic  pair  interactions  and  with 

no  external  magnetic  field,  the  Hamiltonian  H  ,  may  be  taken  to  be  -izza.  .x.  x. , 

n  l  ij  1  j 

where  a.  .  s  0.  The  Curie-Weiss  model  assumes  that  a.  .  =  —  for  all  i  and  j,  that 
ij  ij  n 

is  to  say  that  each  spin  interacts  equally  with  every  other  spin  with  strength 

_1 
2' 


and  takes  P  to  be  symmetric  Bernoulli,  i.e.,  P({-1})=  P((l))  =  y.  Replacing 


P  by  PT(x)  =  P(x/f  ),  we  get 


(2.1) 


dQn(x)  =  z”1  exp[s^/2n]  IldP(Xj), 


where  sn  =  x^  +  . . .  +  xn.  This  model  has  the  advantage  that  the  limiting  free 
spin  per  state  can  be  solved  exactly.  The  existence  of  the  critical  temperature 
and  phase  transition  for  this  model  was  demonstrated  by  Kac  (1968)  .  The 
asymptotic  distribution  for  the  total  magnetism,  Sn,  for  this  model  was  obtained 
by  Simon  and  Griffiths  (1973).  This  is  contained  in  Theorem  2.1  . 

Theorem  2 . 1  (Simon  and  Griffiths).  Let  ,  j  =  1,  . . . ,  n  be  a  triangular  array 

of  random  variables  whose  joint  distribution  is  given  by  (2.1)  and  P  be  symmetric 

3/4 

Bernoulli.  Then  Sn/n  converges  in  distribution  to  a  random  variable  whose 
density  function  is  proportional  to  exp(-y4/l2). 

Theorem  2,1  was  extended  to  the  class  of  probability  measures  L,  which  is 
defined  below,  by  Ellis  and  Newman  (1978b). 

Definition  2,2.  Let  L  be  the  class  of  probability  measures  P  on  R  such  that 

00 

(2.2)  /  exp(x2/2)  dP(x)  <  ®. 

_oo 

Fix  P e  L.  It  can  be  shown  that  condition  (2.2)  guarantees  the  existence 

of  the  moment  generating  function  (m. g.f.),  m(u),  of  P.  Let  h(u)  =  log  m(u)  be 

2 

the  cumulant  generating  function  (c.g.f.)  of  P.  The  function  G(u)  =  u  / 2  -  h(u) 
plays  an  important  role  in  Theorem  2.6  below. 


"■■■I 

--'.--'.I 
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Definition  2.3.  A  real  number  m  is  said  to  be  a  global  minimum  for  G  if 
G(u)  £  G(m)  for  all  u. 

Definition  2,4.  A  global  minimum  m  for  G  is  said  to  be  of  type  r  if 

(2.3)  G(u+m)  -  G(m)  =  c2r  u2r/(2r)  !  +  0(|u|2r+  *)  as  u-H), 

(2r) 

where  c2r  =  G  (m)  is  strictly  positive. 

Definition  2.5.  A  probability  measure  P  is  said  to  be  pure  if  G  has  a  unique 
global  minimum. 

Let  Y  ,  ril,  be  a  sequence  of  random  variables  with  density  function  pr(y) , 

where 

dr  exp  [-c2ry2r/(2r) !  3  if  r  a  2 

(2.4)  Pr(y) = 

N(0,  (1  -  c2) /c2)  if  r  =  1 

and  where  d^  is  the  appropriate  normalizing  constant.  With  these  definitions 
and  notation  we  are  now  in  a  position  to  state  the  generalization  of  Theorem  2.1, 
due  to  Ellis  and  Newman  (1978b) . 

Theorem  2.6.  (Ellis  and  Newman).  Let  P c  L.  Let  P  be  pure,  that  is,  let  m  be  the 

unique  global  minimum  of  type  r  for  G.  Let  ,  j  =  l,  — ,  n,  be  a  triangular 

array  of  random  variables  with  joint  distribution  given  by  (2.1).  Let 

S  =  xjn)  +  . . .  +  X(n) .  Then 
n  1  n 

(2.5)  (Sn-nm)/n1_1/2r  t  Yr, 

where  Y  is  a  random  variable  with  density  function  given  by  (2.4). 
r 

It  is  easily  verified  that  the  symmetric  Bernoulli  measure  is  pure  and 
belongs  to  the  class  L  with  the  corresponding  value  of  r  equal  to  2.  Thus 
Theorem  2 . 6  contains  Theorem  2.1. 
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Note  that  the  moment  generating  function  M(z)  of  the  standard  normal  is 
2 

given  by  exp  (z  / 2).  Thus  we  can  write  (2.1)  as 


(2.6) 


dQn(x)  =  z"1  [M(sn/n)]n  ndP(x^) 


One  might  ask  the  question  whether  it  is  possible  to  obtain  limit  theorems 
of  the  type  (2.5)  when  [M(z)]n  is  replaced  by  the  m.g.f.  4>n(z)  of  a  random  variable 
T  ,  satisfying  some  conditions.  We  answer  this  question  in  the  affirmative  in  the 
next  section. 

3.  Further  Extensions  of  the  Curie-Weiss  Model. 


In  this  section  we  propose  to  extend  Theorem  2.6  by  enlarging  the  class  of 

Hamiltonians  as  well  as  the  class  of  probability  measures  L.  The  large  deviation 

local  limit  theorems  for  arbitrary  sequence  T  ,  n  £  1,  of  random  variables  of  Chaganty 

and  Sethuraman  (1985)  (stated  below)  plays  a  key  role  in  this  extension.  The 

Hamiltonians,  Hn,  in  our  generalized  model  (3.13)  are  taken  to  be  the  cumulant 

generating  functions  of  these  random  variables  Tn< 

Let  (Tn,  nil)  be  a  sequence  of  non-lattice  valued  random  variables  with 

m.g.  f.  's  <(>n(z)  ,  n  >  1,  which  are  analytic  and  non-vanishing  for  z  in  ft  =  { z :  |  Real(z)  |  <  c } , 

with  0  <  c<®.  Let  I  =  (-a,  a)  and  n  =  (z  :  (  Real(z)  j  <  a} ,  where  0  <  a  <  c.  Let 

Si 


(3.1) 


<l>n(z)  =  n  log  ^n ^  *  for  z  e  n 


(3.2) 


Yn(u)  =  sup  [us  -  Vs)]  ,  for  u  e  R. 
I  s  I  <  c 


Let  =  (>|^(s);  sell.  For  uc  An,  we  have  Yn(u)  =  [usn  -  i^n(sn)  ],  where  sn e  I  satisfies 
i|»^(sn)  =  u.  Let  P  be  a  probability  measure  which  satisfies  the  following  condition: 


(3.3) 


/  exp[t^n(x)]  dP(x)  <  oo  for  all  nil, 


Let  h(u)  denote  the  c.g.f.  of  P.  It  is  easy  to  check  that  condition  (3.3)  implies 
that  h(u)  is  finite  for  u  e  R  ,  where 


(3.4) 


Rn  3  <u:  Yn(u)  <  00 }  • 


(3.5) 


Vn(u)  - 


Yn(u)  -  h(u)  for  u  €  Rn 

0°  for  u  4  R  . 

n 


■  '*  / -'y m‘ 


The  function  Vn  plays  the  same  role  as  the  function  G  of  Section  2. 

Definition  3.1.  Let  L*  be  the  class  of  all  probability  measures  P  on  (-c,  c) 
satisfying  condition  (3.3).  We  assume  that  there  exists  L,  p^  > 0  such  that 

Pi 

(3.6)  /  exp[-£V  (u)]  du  =  0(n  x)  , 

*n 

and  V^'s  have  a  unique  global  minimum  at  some  point  mn.  Furthermore  there 
exists  rij  >  0  such  that 

(3.7)  1x1 f  CWu)  ’W3  =  CVn(ran±<S)  '  Vn(mn)]  for  a11  0<<S<rV 
|u|  >  6 

Remark  3.2.  Condition  (3.7)  is  used  mainly  in  inequality  (3.27)  of  Lemma  3.13. 

An  easily  verifiable  sufficient  condition  for  (3.7)  is 

(3.8)  V'(u)  >  0  for  u  >  m  and  V'(u)  <  0  for  u  <  m  . 

n  n  n  n 

In  all  the  examples  of  section  4  we  will  be  verifying  (3.8)  instead  of  (3.7). 

Remark  3.3.  Suppose  that  =  (-»,«).  If  Yn(u)/|u|  converges  to  ®  as  |u|  ■*■<*>, 
then  condition  (3.3)  implies  (3.7)  as  seen  below: 

(3-9)  exp[-Vn(u)]  =  exp[--yn(u)  +  h(u)] 

=  expC-y  (u)][  /  exp[ux]  dP(x)  +  /  expLux]  dP(x)] 

|  x  |  <  A  I  x  f  >  A 

^  exp[-y  (u)  +  uA]  +  /  expC^  (x)]  dP(x) 

"  l*[>A 

Y  (n) 

<  exp[-|u|  (———  -  A) ]  +  /  exp(4in(x)]  dP(x). 

1  1  !  x  |  >  A 

The  right  hand  side  can  be  made  close  to  zero  first  by  choosing  A  and  then  letting 
|u|  -*•«•.  This  shows  that  Vn(u)  ■*■<*>  as  |u|  -*■«.  Since  mn  is  the  unique  global  mini¬ 
mum  of  V  ,  this  also  shows  that  condition  (3.7)  holds. 
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Let  m  e  A  .  Then  there  is  a  t  in  I  such  that  ^'(t)  =  m  .  For  tel, 
n  n  n  n'  n  ’ 

define 

(3.10)  Gn(t)  =  ^n(Tn)  +  itmn  -  i|>n(-rn  +  it)  . 

The  following  theorem,  which  provides  an  asymptotic  expansion  for  the  density- 
function  kn  of  Tn/n  in  terms  of  the  large  deviation  rate  yn ,  is  due  to  Chaganty 
and  Sethuraman  (1985) . 


Theorem  3.4.  Assume  the  following  conditions  for  T^ : 


(A)  There  exists  8  >  0  such  that  |  ip  (z)  |  <  8  for  z  e  ft  and  n  2  l, 

n  cL 

(B)  There  exists  a  >  0  such  that  <J>^(t)  £  ct  for  t  e  I  and  n  2:  1. 

(C)  There  exists  n  >  0  such  that  for  any  0  <  5  <  n, 

inf  Real  (G  (t))  =  min[Real(G  (6)),  Real(G  (-6))],  for  nil. 
|t|>6  n  n  n 


(D)  There  exists  p  >  0  such  that 

oo 

sup  /  Un(T  +  it) /4>n (t)  1 1/n  dt  =  0(np). 

T  £  I 

Then  1/2 


(3.11) 


k  (m  )  = 
nv  n' 


2ir*nCTn3 


exp(-nyn(mn))[l  +  O(-) ]. 


Remark  3.5,  When  T^  is  the  sum  of  n  i.i.d.  random  variables,  condition  (C)  is 
automatically  satisfied  and  conditions  (A),  (B) ,  (D)  are  easy  to  verify,  since 
they  do  not  depend  on  n . 

Remark  3.6.  Suppose  that  m  e  nA  .  Then  there  exists  Cel  such  that  iJi"(E  )  =  m, 
-  vv  n  n  n  n 

for  n  2  1.  In  this  case  one  can  verify  that  vV  f E  )/  /ii>"  (t_ )  =  [l+0(lm  -ml)] 

n  n  n  n  n 

and  thus  we  can  rewrite  (3.11)  as 


The  cumulant  generating  function  of  P  is  given  by 


(4.6)  h(u)  =  log(3/2)  -  3  log|u|  +  |u|  +  log[|u|  (l  +  e““'ul)  -  (l-e"2'u')3,  ue  R. 


Therefore  u  e  R, 


(4.7)  V(u)  =  y(u)  -  h(u) 

=  C- 1  +  ^1+u2  ]  +  log|u|  +  log[-l+  A+u2  ]  -  |u| 

-  log[|u|(l  +  e'2lul)  -  (l-e"2lui)]  +  log(4/3) 

Part  (I) ,  We  just  need  to  verify  (4.1).  Now 

(4.8)  V  ' (u)  =  A+(  1/u2T  +  (2/u)  -  u/(ucoth(u)  -  1) 

>  1  +  (2/u)  -  u/(ucoth(u)  -  1) 

-  (u2  ♦  u)  (coth(u)  -  1)  +  ucoth(u)  -  2  >  Q 
(ucoth(u)  -  1) 

for  u  2  2,  since  ucoth(u)  -  1  >  0  and  coth(u)  >  1.  Again 


(4.9) 


y.(u)  =  il-Uii  for  u  >  0 

u[ucoth(u)  -  1] 


By  differentiation  one  can  verify  that 


(4.10) 


2  4 

ucoth(u)  -  1  >  ^ for  0  <  u  <  7.5 


(4.11) 


2  4 

u  u 


3  45  2  +  /l  +  u2 


for  0  <  u  <  2.19. 


These  two  inequalities  show  that  V'  (u)  >  0  in  the  region  0  <  u  <  2.  Since 
V”  (0)  =  y"  (0)  -  h"  (0)  =  j  -  ,  V  has  a  unique  global  minimum  of  type  1. 


*  -  ■*,  ■’>  ;/^v.  .  -v  •*.  ■*. 


written  down.  Part  (I)  will  show  that  V  has  a  unique  global  minimum  of  type  r  at 
mn:0.  Part  (II)  will  verify  the  four  conditions  of  Theorem  3.4  or  Theorem  3.16 
depending  on  whether  Tn  is  non-lattice  or  lattice.  Part  (III)  will  verify  condi¬ 
tions  (3.3),  (3.6)  and  (3.8)  that  insure  that  P  e  L*,  in  the  case  of  continuous 
P  and  conditions  (3.3),  (3.8)  and  (3.38)  that  insure  that  P  e  L^,  in  the  case  of 
discrete  P.  These  three  parts  will  impiy  that  the  random  variable  S  ,  after  proper 
normalization  converges  in  distribution  to  Y*,  for  appropriate  value  of  r. 

The  verifications  of  parts  (II)  and  (III)  are  usually  very  easy.  The  veri¬ 
fication  of  Part  (I)  presents  more  difficulties.  However,  in  all  the  examples 
considered.  Vis  symmetric  around  the  origin  and  hence  Part  (I)  is  verified  once 
we  show  that 

(4.1)  V'(u)  >  0  for  u  >  0. 

Note  that  (4.1)  also  verifies  condition  (3.8).  The  verification  of  (4.1) 
is  routine  in  all  the  examples,  but  sometimes  long,  and  therefore  only  the  main 
steps  are  presented.  One  should  also  note  that  in  part  (II),  condition  (C)  or 
(C1 )  is  automatically  satisfied  by  Remark  3.5. 

Example  4.1.  Let  F  have  p.d.  f.  (1/2)  exp(-|x|),  -«<x<®,  and  P  have  p.d.f. 

(3/4) (1  -  x  )  for  |x|  <1.  Thus  the  joint  distribution  given  in  (3.13)  becomes 

(4.2)  d^Oc)  =  z~1(3/4)n(l  -  sj/n2)‘n  JI(l-x?)  Mx. . 

In  this  case 

(4.3)  <|>n(z)  =  U  •  z^)  and 

(4.4)  4>(z)  =  -  log  (1  -  z^ )  for  |  Re  al  ( z )  j  <  1, 

(4.5)  y(u)  =  log  2  +  [-1  +  vl  +  u2  3-2  log  |u|  +  log  [-1  +  A+u^],  for  u  e  R. 


where 


<(y„)  ■  (i.Vr)/2  kn("n*n'1/2\)  exp[n(h(mn*n-1/2ryn).vn(inn))3. 


Imitating  the  proofs  of  Lemmas  3.9  thru  3.13,  one  can  show  the  following: 


1  *- 


(i)  - gn(y  )  "*■  g(y)  as  n  +*,  uniformly  on  bounded  intervals  of  y. 

K\ 


(ii) 


■I  l^ri^n1  "  /  g(y)dy  aS  n 

s  n  -°° 


(Ui)  ,  ,?  l/4r*W  *  0 


as  n  -*•  «. 


Iynl>n 


The  above  three  steps  (i) ,  (ii) ,  (iii)  complete  the  proof  of  Lemma  3.19. 


4.  Applications. 

In  this  section  we  illustrate  the  main  theorems  of  section  3  with  four 

applications  and  demonstrate  limit  theorems  in  quite  complicated  situations  of 

dependent  variables.  The  model  (3.13)  for  the  joint  distribution  of  (X^  , . . .  ,X^n^) 

is  completely  specified  if  Tn  and  P  that  arise  in  it  are  specified.  To  simplify 

matters,  in  all  the  examples  of  this  section  we  let  T^  be  the  sum  of  n  i.i.d. 

random  variables  with  common  d.f.F.  The  four  examples  below  contain  all  occurences 

of  lattice  and  non-lattice  T  ,  and  continuous  and  discrete  P.  The  limit  distri- 

n 

bution  of  the  normalized  sum  S  =  X^n^  +  ...  +  X^  is  normal  (r=  1,  in  the  notation 

n  1  n 

of  Theorems  3.7  and  3.18)  in  example  4.1  and  is  non-normal  (r=2)  in  examples 
4.2,  4.3  and  4.4.  The  results  of  Ellis  and  Newman  (1978b)  show  that  limit  distri¬ 
butions  with  every  possible  value  of  r > 2  can  also  arise  in  suitable  models. 

In  all  the  examples  below  we  will  specify  F  and  P  and  write  down  the  joint 
distribution  Q^.  The  exact  expression  for  <f>n,  ij»n  =  <p,  Yn  =  Y»  h  and  vn  =  V  will  be 


If  T  satisfies  the  conditions  of  Theorem  3.16,  then 
n  ’ 

(3.40)  (S  -  nx  )/n1_1/2r  +  Y*  , 

n  n  '  r  ’ 

where  Y*  and  r  are  as  defined  in  Theorem  3.7. 
r  n 

The  proof  of  the  above  theorem  parallels  the  proof  of  Theorem  3.7.  We 
therefore  outline  briefly  the  modifications  that  need  to  be  done.  Note  that 
dQ^  can  be  written  as 

(3.41)  d^U)  =  l  IUMn  y(x.)  f*(y) , 

where  f*(y)  =  z”1  kn(mn  +  n”1//2ry)  exp[n  h(m^  +  n_1^2ry)  3  is  a  probability  mass 

function  of  a  lattice  valued  distribution  with  span  h'  =  h  /n^_1^2r,  and  dM  (x.) 

v  n  n  ’  n,y  y 

is  as  defined  in  (3.30).  We  introduce  discrete  random  variables  W*  with  p.m.f.  f*. 
It  suffices  to  show  that  W*  converges  weakly  to  a  continuous  random  variable  W 
with  probability  density  function  f,  defined  in  (3.35).  The  rest  of  the  proof 
is  identical  to  the  proof  of  Theorem  3.7.  Note  that  the  span,  h,i,  of  W*  converges 
to  zero.  By  a  theorem  of  Okamoto  (1959),  the  sequence  of  random  variables  W* 
will  converge  in  distribution  to  W,  once  we  prove  the  following: 

Lemma  3. 19.  For  yeR,  define  yn  =  h^[y/h^3.  Let  the  probability  mass  function 
fn  and  the  probability  density  function  f  be  as  defined  above.  Then 

(3.42)  — ~  f*(yJ  f(y)  as  n  +*, 

K\ 

uniformly  on  bounded  intervals  of  y. 

Proof  (outline).  Note  that  f(y)  =  g(y)//g(y)dy,  where  g(y)  is  as  defined  in 
(3.15).  We  first  write 
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(C')  There  exists  n  >  0  such  that  for  any  0  <  5  <  n, 


inf  Real(G„(t))  =  min[Real(G„  (<5) ) ,  Real(G  (-6))]  for  nil, 

«  *  |t|  s  ir/|hn|  n  n  n 

where  Gn(t)  is  defined  by  (3.10). 

(D')  There  exists  p  >  0  such  that  Ih^j-1  =  0(np). 


As  before  for  a  probability  measure  P  on  R,  define  V^fu)  as  in  (3.5). 
The  class  of  probability  measures  that  are  of  interest  is  defined  below. 

Definition  3. 17.  Let  L*  be  the  class  of  probability  measures  P  satisfying 
conditions  (3.3),  (3.7)  and  (3. 38) (defined  below) . 


(3.38)  £  exp[-£V  (u)]  =  (Kn^1)  for  some  £,  p..  >0. 

ueR  n  1 

n 

Note  that  (3.38)  is  the  appropriate  replacement  of  (3.6)  for  the  lattice  valued 
case. 

For  Hamiltonians  which  are  functions  of  the  moment  generating  functions  of 
lattice  valued  random  variables  we  have  the  following  theorem  almost  identical  to 
Theorem  3.7. 


Theorem  3.18.  Let  P  e  LJ  . 
random  variables  satisfying 


Let  X.n  ,  j  =  l,  ...»  n,  be  a  triangular  array  of 
|  <  c  and  having  a  joint  distribution  given  by 


(3.39) 


dVx)  =  z'n  Vsn/n)  ndp(xj)« 


where  *  is  the  m.g.f.  of  the  lattice  valued  random  variables  T  .  Let 
n  n 

Sn  *  x|n)  ♦...•*•  x<n) .  Let  Vn  have  a  unique  global  minimum  of  type  r  at  the 
point  ®n  e  An.  Let  converge  to  a  point  m  belonging  to  the  interior  of  nAn> 


The  unconditional  limiting  distribution  of  (Sn  -  nTn)/n*~*^r  is  just  the  mixture 
of  the  limiting  conditional  distribution  and  f(y),  by  Theorem  3.15  of  Sethuraman. 
This  completes  the  proof  of  Theorem  3.7.  ]| 

Remark  3.14.  When  Tr  is  the  sum  of  independent,  normally  distributed  random 

2 

varibles  with  mean  zero  and  variance  one,  4>n(sn/n)  becomes  exp[sn/2n]  and  the 
class  of  probability  measures  L*  reduces  to  the  class  L.  Thus  Theorem  3.7 
generalizes  Theorem  2.6  to  a  larger  class  of  Hamiltonians  and  Probability  measures. 

We  now  state  the  theorem  of  Sethuraman  (1961)  which  was  crucially  used  to 
obtain  the  limiting  marginal  distribution  of  (Sn  -  nTn)/n*-1^r  in  the  proof  of 
Theorem  3.7. 

Theorem  3.15  (Sethuraman).  Let  Ar  be  a  sequence  of  probability  measures  on  V*W, 
where  V  and  W  are  topological  spaces.  Let  un  be  the  marginal  probability  measure 
of  on  V  and  vn(v,-)  be  the  conditional  probability  measure  on  W.  Suppose  that 
un  converges  to  a  probability  measure  u  for  every  measurable  set  in  V  and  for 
almost  all  v  with  respect  to  u,  v  (v,»)  converges  weakly  to  v(vf.).  Then  Ar 
converges  weakly  to  A,  where 

(3.36)  A(A*B)  =  /  v(v,B)  du(v) 

A 

for  every  measurable  rectangular  set  A*B. 

We  now  turn  our  attention  to  the  case  where  T  ,  n>  1,  are  lattice  valued 
random  variables  with  spans  hn,  nil.  The  following  theorem,  which  is  analogous 
to  Theorem  3.4,  was  proved  by  Chaganty  and  Sethuraman  (1985). 

Theorem  3.16.  Let  T  ,  n  21,  be  a  sequence  of  lattice  valued  random  variables  with 
spans  hn,  nil.  Let  mn  belong  to  the  range  of  Tn/n.  Assume  that  conditions  (A), 
(B)  of  Theorem  3.4  hold  and  replace  conditions  (C) ,  (D)  by  the  following: 
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We  first  note  that 


(3.32)  log  Em  expCt(Sn  -  nt^/n1  '  1/2r) 
n,y 


=  n 


n 


l-l/2r 


+  h 


vm-"  +  m  (y) 

l-l/2r  n,rw' 


h("n,r(),)) 


*  n 


CT757  *  h'tmn,rM)  *  h''  ^ttf  *  "l"'1) 


=  h"  (m  )  ty  + 


2  n 


17171 


+  o(l)  , 


since  tn  =  h'(mn).  Thus 


(3.33)  log  Em  exp[t(Sn  -  n-tn)/n1_1^2r] 
n,y 

h"(m)ty  if  r  2:  2 
h"(m)ty  +  h— if  r*  1. 


This  shows  that  the  limiting  distribution  of  (S^-n-r^/n1  ^2r  given  Wr  =  y 
is  degenerate  at  h"(m)y  if  r  >  1  and  N(h"(m)yf  h"(m))  if  r  =  1.  Next  we  note 
that 


(3.34) 


fn(y)  =  z'1  n'1/2r  kn(mn  rCy) )  exp(n  h(mn  r(y))) 


gnCy) /  /  ^(y)  dy 


where  gn(y)  is  as  defined  in  (3.16).  By  lennnas  3.9,  3.11  and  3.13  it  follows 
that 

(3.35)  f  (y)  -  f(y)  =  ■  as  n  -*■  ®  , 

/g(y)dy 

where  g(y)  =  exp[-y  C2r/(2r)!].  Thus  the  limiting  distribution  of  W^  is  f(y) . 
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Hence 

lwL^g"w  dy 


=  OCn^)  exp[-(n-£)  [■ 


c2r,n  J_ 
(2r) !  v'n 


+ 


E  n-(2r*l)/4r]] 


which  goes  to  zero  since  | |  <  K  for  all  n.  The  proof  of  Lemma  3.13  is  now 
complete. 

Proof  of  Theorem  3.7.  We  first  express  dC^  defined  in  (3.13)  as  follows: 

(3.28)  dQn(x)  =  z^1  <frn(sn/n)  ndP(Xj) 

■  z"1  /ex p(ysn)  kn(y)  dy  HdP(x^) . 

_  1  /  2  T 

Substituting  mn  r(y)  *  mn  +  n  y,  we  get 

(3.29)  dQn(x)  =  z^1  n"1/2r  /ex p(mn>r(y)  sn)  knCmn>r(y))  dy  HdP(x.) 

■  z"1  n‘1/2r  /n  exp(xj  mn >r(y)  -  h(mn>r(y))  dP(Xj) 

*  kn(mn,r(y))  exp(nh(n^  (y))  dy 
=  /ndMn  y(Xj)  fn(y)  dy,  where 

(3.30)  dMn  y(Xj)  =  exp(Xj  mn^r(y)  -  h(mn  p(y)))  dP(x^)  and 

(3.31)  fn(y)  ■  z^1  n"1/2r  kn(mnjI.(y))  exp(n  h(n»n  r(y))) . 


Since  / do  (x)  =  1  and  /dM  (x.)  =  1  for  each  y  and  j,  we  have  / f  (y)dy  =  1. 
ti  n  ,y  j  n 

we  can  introduce  random  variables  with  probability  density  function  f^(y) 
the  representation  (3.29)  of  dQ^x)  shows  that  given  wn  =  y»  ,  j  =  1,  ..., 
i.i.d.  with  common  distribuiton  M  (x) . 

“  *y 


We  now  proceed  to  obtain  the  limiting  distribution  of  (Sn  -  mn)/n 


l-l/2r 


dM  (x) . 

n.y 


Thus 
and 
n  are 

under 


Substituting  z  =  n 


l/2r 

y,  we  get 

<  S2irte  )  n'(1'1/r)/2  ,  ,/  -i/4r|expC‘n(VVz)  "  Vmn))]1 

ti  n  J  2  J 

•  |  exp(n  Yn(n>n  +  z))  kn(mn  +  z)|dz 

<  0(np+U+1/r)/2)^max^_1/4r  exp[-(n-£)  (Vn(mn  +  z)  -  Vn(mn))] 

•  /exp(-£(Vn(mn+ z)  -  Vn(mn)))dz 

The  last  inequality  follows  from  Lemma  3.12.  This  together  with  condition 
yields 


^  °(nC1)  |zf>n~1/4r  expC*(n"e)(Vn(mn  + z)  ‘  Vm, 
*  0(nq)  exp[- (n-£) L^] , 

where 


(3.27) 


iLl/4r*n(rt  ■'>' 


q  =  Pl  +  p  +  (1  +  l/r)/2, 
and 


z 


±  n 


•l/4r 


|^n-V4rCVn(“n*z)  -  W5' 

by  condition  (3.7).  Therefore, 


This  minimum  is  attained  at 


Ln  =  minC(Vn< 


+  n" 1/4r) 


W>* 


<VMn 


n'1/4r)  -  W” 


c2r,n 

'  127)7 


1  + 


K 


n 


n 


(3.6) 


))] 


-  (2r  +  l)/4r 
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(3.24)  sup  |exp(n  yn(mn  +  y))  kn(mn  +  y)|  =  0(np+1)  as  n  ■*  “  . 

Proof.  An  application  of  the  inversion  formula  yields  (see  (2.12)  of  Chaganty 
and  Sethuraman  (1985)), 

|exp[n((mn  + y)s  -  tf>n(s))3  kn(mn  +  y)| 

00 

=  \y-  !  exp[n(ij)  (s  +  it)  -  ip  (s)  -  it(m  +y))3  dt| 


5  Tz  I  Un(s  +  it)  /  $ _ (s)  1 1/n  dt  . 


2ir  J  ,rn 

-oo 


Taking  supremum  with  respect  to  s  e  I  and  using  condition  (D)  of  Theorem  3.4  we 

get  . 

sup  |  exp(ny  (m  +  y) )  k  (m  +  y)  |  =  0(np  )  .  II 


Lemma  3.13,  Suppose  that  Vn  has  a  unique  global  minimum  at  the  point  mn  e  An 


and  let  gR  be  as  defined  in  (3.16).  Then 


(3.25) 


.  ,/  1/4r  £,(>0  dy  +  0  as  n  -  » 

|y|>n 


Proof.  Let  m^  r(y)  =  mn  +  n  *^ry.  By  (3.16)  we  have 

r2,  *-(.)] 1/2 

r  _  j..  _  n  r 


(3.26) 


\y\>n 


/  , . .  k  (m  (y)) 
|  £  nl/4r  n,rv/  J 


exp[n(h(mn^r(y))  +  Vn(mn))]  dy 


2»  *n  (  V 


Jl/4r  kn(mn,rty^ 


exp  [C-n(Vn(mn>r(y))  -  Vn(mn))  3  +  nyn(mn  rCy))3  dy 
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Proof.  Note  that  n  *y2ry  converges  to  zero  uniformly  in  y  for  jy|  <  n^4r. 

Since  m  is  an  interior  point  of  nA  there  exists  N.  (independent  of  y)  such  that 

n 

-  l/2r 

mn  (y)  =  (mn  +  n  y)  e  Afl  for  n  £  N^.  Applying  Theorem  3.4  for  n>N3,  we  get 

1/2 

2 it  i|t"  (E  ) 

(3<22)  ,  l/4r  *n(y)dy*  - n -  ,  ./  l/4r  exp[n(h(i^  (y))  +  v„(\)) 3 

|y|sn  L  J  |y|-n 

*  kn(ran,r(y))  dy 

=  .  J  l/4r  exPC'n(VnCmn,r(y))  '  W}1 
I  y  I  — n  ' 

•  [1  +  0(!mn  r(y)  -  m|)  + 0(i)]dy 

00 

■  /  An(y)  dy  , 


where 


*n(y)  «  I ( |y|  s  n  1/4r)  exp[-n(Vn(mn>r(y))  -  vn(mn))] 

•  [1  +  0(|mn^r(y)-roj)  + 0(^-)]  , 


and  I(*)  is  the  indicator  function.  It  follows  from  Lemma  3.10  that  |An(y)|  is 
bounded  by  an  integrable  function.  We  can  now  conclude  from  Lemma  3.9  and 
Lebesgue  dominated  convergence  theorem  that 


(3.23) 


/  *n(y)  dy  -*■  J  g(y)  dy  as  n 


The  proof  Lemma  3.11  is  now  complete. |j 

The  next  Lemma  3.12  is  needed  in  the  proof  of  Lemma  3.13. 


Lemma  3.12.  Let  T  ,  n  2  1  be  a  sequence  of  random  variables  satisfying  the 
conditions  of  Theorem  3.4.  Then 


(3.20) 


nCvn(mn  +  n‘1/2ry)  -  V^n^)]  >  y2r  c2r/2(2r)! 


for  all  n  s  N,  and  |y|  <  n^4r. 


Proof.  Let  0  <  e  <  c2r/2.  Since  c2r  n  converges  to  c2r  we  can  find  such 

that  c_  >  c-  /2  +  €  for  all  n^N,.  Recall  that  y  (m  )  *  m  t  -  \p  (t  ) ,  where  t  is 
2r,n  2r  1  'n  n  n  n  rnv  n  *  n 

such  that  ^ ( tr )  =®n>  It  is  easy  to  verify  that  Y^(mn)  =  Tn  and  y”  (m^)  =  Ci|>"  (Tn)  • 

Also,  Y^2r+1^(ran)  the  ratio  of  a  simple  function  of  (2r+l)  derivatives  of 
i|>n  at  xn  to  Cip '*  (^)  ]2r  .  Conditions  (A)  and  (B)  of  Theorem  3.4  imply  that  all 
these  derivatives  are  bounded  uniformly  in  n  and  that  <j i"(x  )  £  a  >  0  (see  (2.6) 
of  Chaganty  and  Sethuraman  (1985)).  Hence  bounded  in 

n  and  consequently  \f^T+^  (mn)  =  Y^2r+^  (®n)  -  h^2r+1^ (mn)  is  also  uniformly 
bounded  in  n.  Therefore 

W“>  -  W  *  “2r  W(2r)!  *  V2"1 

as  u  +  0,  where  (K^j  s  K  <  <*>  for  all  n.  Thus 

n[V  (m  +  n  ^2ry)  -  V  (m  ) ] 
n  n  1  nv  n' 

2r  v  2r+l .  l/2r 

=  y  c2r,n/(2r)!  +  Kn  y  /n 

^  y2rc2r/2(2r)  !  +  y2r[e/ (2r) !  -  K  y/n^2*] 

>  y2rc2r/2(2r)l 

if  (y|  <  n^4r  and  n  £  N  =  max{N^,  (K(2r)  !/e)4r) .  This  completes  the  proof  of 
the  Lemma  3.10.  || 

Lemma  3.11.  Let  g  and  g^  be  as  defined  in  (3.15)  and  (3.16).  Then  under  the 
hypothesis  of  Lemma  3.9  we  have 


(3.21) 


,  .  /  1/4r  «n(y)  dy  ->  /  g(y)  dy  as  n 
| y |  5n  -® 
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For  y  c  R ,  let 


(3. IS) 


g(y)  =  wC-y  c2r  /  (2r)  !  ]  and 


-l/2r  -l/2r 

(3.16)  ^(y)  =  I- - „ - J  kn(mn  +  n  1/2  y)  exp[n(h(mn  +  n  y)  +  Vn(mn))3 

where  Cn’s  are  defined  as  in  Remark  3.6.  The  functions  g^s  arise  in  the  proof 
of  Theorem  3.7.  Lemma  3.9  shows  that  g^(y)  converges  to  g(y)  as  n->-“  for  each 
y.  The  next  four  lemmas,  Lemma  3.10  -  3.13,  show  that 


(3.17) 


/  gjjCy)  dy  /  g(y)  dy  as  n 


Lemma  3.9.  Suppose  that  has  a  unique  global  minimum  of  type  r  at  the  point 

mn  e  A^.  Let  mn  converge  to  m,  where  m  is  an  interior  point  of  nAn.  Suppose  that 
c  2rl 

\r  J  (m  )  =  c,,  converges  to  c~  as  n -*•«>.  Then 

n  v  nJ  2r,n  °  2r 


(3.18) 


gn(y)  ■*  g (y)  as  n  -*■  ®  . 


Proof.  Fix  y  €  R.  Let  mn  (y)  =  m  +  n'1/2ry-  Then  m  (y)  converges  to 
— -  n  j  r  n  ^  ^ 

m  and  mn  r(y)  e  An  for  sufficiently  large  n.  Applying  Theorem  3.4  together  with 
Remark  3.6,  with  mn  replaced  by  mn  r(y)  we  get 

(3.19)  gn(y)  =  exp[-nyn(mn^r(y))  +  n(h(mn>r(y))  + 

•  Cl  +  0(|mn  (y)  -  ml)  +  0(^)3 

=  exp[-n(Vn(mn  r(y))  -  Vn(mn))]  [1  +  0(imn  r(y)  -  m[)  +  0(^)3 
=  exp[-y2rc2r  n/(2r) !  +  n  o( | y| 2r/n) 3  [1  +  0(  |mn^r(y)  -  m| )  +  O(i-)  3 
+  g(y)  as  n  -*•  »  .  || 

Lemma  3. 10.  Suppose  that  the  Vn’s  have  a  unique  global  minimum  of  type  r  at  the 

point  m  e  A„ .  Then  there  exists  N  such  that 
r  n  n 
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(3.12) 


k  (m  )  = 
n  n 


1/2 


2ir 


-ny  (m  )  , 

!  [1  +  0(jmn  -  m| )  +  0 (“} 3 . 


For  each  integer  rsl,  let  Y*  be  a  random  variable  with  probability  density 
function  given  by  dr  exp[-C2T  y2r/[h"  (m)  32r  (2r) !  ]  if  rs2  and 
N(0,  h"  (m)  [h"  (m)  +  C23/C23  if  r=  1,  where  C2r  is  the  constant  that  appears  in 
Theorem  3.7  below  and  df  is  the  normalizing  factor.  With  these  assumptions 
and  notation,  we  are  in  a  position  to  state  the  main  theorem  of  this  section. 


Theorem  5.7.  Let  ,  j  =  1,  . . . ,  n  be  a  triangular  array  of  random  variables 


satisfying  | Xj 


(3.13) 


(n). 


<  c  and  having  a  joint  distribution  given  by 
-1 


dQ^x) 


=  z 


n 


wn) 


ndP(x.), 


where  <Pn  is  the  m.g.f.  of  Tn  and  P  e  L*.  Assume  that  Vn,  defined  in  (3.5),  has 

(2r) 

a  unique  global  minimum  of  type  r  at  mn  e  Afl.  Let  n^-nn  and  \r  (ran)  '►^r  as 

n  •>  ”,  where  m  is  an  interior  point  of  nA  .  Let  +  ...  +  .  If 

’  r  n  n  1  n 

T  satisfies  the  conditions  of  Theorem  3.4,  then 
n 

(3.14)  (Sn-nxn) /n1_1/2r  £  Y*  , 


where  ip '  ( t  )  =  m  and  Yi  is  as  defined  above. 
nv  n  n  r 

The  proof  of  the  above  theorem  is  postponed  until  the  end  of  Lemma  3.13. 


Remark  3.8.  The  distribution  function  (^(x)  is  well  defined  because 

zn  =  /exp  Cn  ipnCsn/n)  3  ndP(x^) 
c 

£  (  /exp  [  ip  (x)  ]  dP(x)]n<°°, 

-c  n 

wherein  we  have  used  condition  (3.3)  and  the  fact  that  ip  is  a  convex  function. 
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Part  (II).  It  is  easy  to  verify  that  the  random  variables  T  ,  ns  1,  satisfies 
conditions  (A),  (B)  of  Theorem  3.4.  Condition  (D)  also  holds  because 


(4* 12) 


00 

sup  /  j  4(x  +  it)/ 4> (x )  |  dt 
|  t  |  <  a  <  1  -® 


< 


SUP 
t  |  <  a 


2  +  2  / 
1 


/  2(1 -t2) 

0  [4t2  t2  +  (1  -  t2+  t2)]^ 

7?  dt  <  “• 


dt 


Part  (III).  To  show  P  e  L*,  we  only  need  to  verify  conditions  (3.3)  and  (3.6). 
Condition  (3.3)  trivially  holds  since  exp[<J>(u)]  is  bounded  on  (-1,  1),  the  support 
of  P.  Condition  (3.6)  on  the  integrability  follows  from  the  fact  V(u)  *v  log  |u| 
as  |  u|  -+  ®. 

We  can  therefore  conclude  in  this  example  that 
(4.13)  Sn/^n  +  N(o,  1/3). 


Example  4.2.  Let  F  be  the  distribution  function  of  the  sum  of  two  independent 
and  identically  distributed  uniform  random  variables  on  the  interval  (-b,b)  with 
b  =  /3//2  .  Let  P  be  standard  normal  probability  measure.  The  joint  distribution 


is  given  by 
(4.14)  dQ^x) 


-1„  ,.„/2rn  sin«bV")  ,2n  rl 
(2,)  C - (ST) -  3  6XP  C’T 

v  n 


ix.  ]  Hdx.. 
i  j 


In  tnis  example 

(4. 15)  $n(z)  =  [sinh(bz)/(bz)]2n 

(4.16)  \p(z)  =  2[log(sinh(bz))  -  log(bz)]  for  |Real  (z)  |  <  ®  . 
Also 

(4.17)  y(u)  =  sup  [us  -  4i(s)  ]  for  |uj<2b. 
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(4.18) 


V(u)  =  Y(u)  -  h(u)  =  Y(u)  -  (u/2) 

=  usx  -  -  (u2/2)  ,  for  |u|  <  2b, 


where  \|/(s^)  =  u. 


Part  (I).  Differentiating  (4.18)  with  respect  to  u  we  get 


(4.19)  V'(u)  -  CSl  -  ♦'(Sj)] 

2 

=  [s^  -  2bs1  cothfbsp  +  2]  /  s^  ,  where  ^ ’  (Sj)  =  u  . 

2 

Taking  successive  derivatives  one  can  show  that  t  -  3(tcoth(t)  -  1)  >  0  for  t  >  0. 
Letting  t  =  bSpand  noting  that  b  =  /l //I  ,  we  get 


(4.20) 


Sj^  -  2bs1  coth(bs^)  ♦  2  >  0  for  >  0. 


Thus  V'(u)  >  0  for  u  >  0.  With  our  choice  of  b  *  /3  / /I  one  can  verify  that 
f41 

V"  (0)  =  0  and  Vv  J  (0)  =  3/5  >  0.  Therefore  V  has  a  unique  global  minimum  of 
order  2  at  the  origin. 


Part  (II).  We  need  to  verify  that  Tn  satisfies  conditions  (A)  thru  (D)  of 
Theorem  3.4.  Note  that  ij;"(0)  =  1.  Hence,  there  exists  a  >  0  such  that 
<1>"(t)  >  j=a  for  |  x  |  <  a.  This  verifies  condition  (B) .  It  is  easy  to  check 
that  condition  (A)  also  holds  for  this  choice  of  a.  Consider 


sup 
t  |  <  a 


(sinh  (bt))2  +  (sin(bt))2 


.2 ,  2  2. 
b  (t  + 1  ) 


C  ,  -00  <  t  <  “  , 


for  some  constant  c.  Therefore 
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(4.21)  /  |  4>(t  +  it)/<j>(t)  J  dt  s  /  ldt+  /  dt  < 


t|  SI  1 1 1  >  1 


This  verifies  condition  (D). 


Part  (III).  We  only  need  to  verify  that  conditions  (3.3)  and  (3.6).  Let 

e  >  0  be  given.  There  exists  a  5  >  0  such  that  | x|  <  6  implies  [sinh(bx)/(bx) ]  <  (1+e).  » 

Therefore 

OO  QO 

(4.22)  /  expO(x)]  dP(x)  =  (2w)  ^2  /  [sinh(bx)/(bx) ]  exp(-x2/2)  dx 

-CO  -  00  ~ 

I 

<  (1  +  e)  +  (2it) ~ 1/2  |x|/fi  Csinh(bx)/(bx)]2  exp(-x2/2)  dx 
s  (1  +  e)  +  (b6)~2C2ir)"1/2  .  J  ^inh(bx)]2  exp(-x2/2)  dx, 

|  X  |  >  6 

t 

which  is  finite.  Hence  condition  (3.3)  is  satisfied.  The  other  condition  (3.6) 

also  holds  since  the  integral  is  over  a  finite  range  and  V(u)  £  -h(u). 

Thus  the  random  variables  T  's  and  P  satisfy  the  conditions  of  Theorem  3.7.  In  T 

n  . 


this  example 


(4.23) 


V"3/4  4  *1  • 


where  the  p.d.f.  of  Y|  is  given  by  exp(-y  /40) 


-00  <  V  <  oo 


Example  4.3.  Let  F  be  as  defined  in  Example  4.2.  Let  P  be  symmetric  Bernoulli, 

i.e.,  P({-1})  =  P( { 1>)  =  ty 2  .  The  joint  distribution  is  given  by 

j—  -r  2n 

,  n  sinh(bs  /n) 

(4.24)  dO  (x)  =  z'1  — - 2 - 

^  -  n  (/lbs) 


where  x^  =  +1  for  all  1  £  j  £  n,  and  b  =  /3  /  /?  .  Recall  that 

(4.25)  4>(z)  =  2  log(sinh(bz))  -  2  log(bz)  for  |Real(z)|  <  «, 


(4.26)  y(u)  =  us1  -  *(Sl)  if  |u|  <  2b, 


<  00 


where  ^'(sp  =  u.  The  cumulant  generating  function  of  P  is  given  by 

(4.27)  h(u)  =  log  cosh(u)  ,  -<*>  <  u  <  «>  . 

The  function  V(u)  =  y(u)  -  h(u)  is  finite  if  u  <  2b.  Now  for  u  >  0, 

(4.28)  V ' (u)  =  sx  -  h •  (i|» •  (Sj) ) 

*  Sj  -  tanh(>p'(s1)) 

>  -  ♦'(Sj)  >  0  , 

as  shown  in  Example  4.2.  Since  the  first  three  derivatives  of  V  at  the  origin  are 

(4) 

equal  to  zero  and  V  (0)  =  13/5,  the  point  zero  is  a  minimum  of  order  2  for  V. 
This  completes  the  verification  of  Part  (I).  We  have  already  checked  Part  (II) 
in  Example  4.2.  One  can  easily  show  that  the  probability  measure  P  belongs  to 
the  class  L*,  completing  verification  of  Part  (III).  Thus  by  the  conclusion  of 
Theorem  3.7  we  get 

(4.29)  V1^4  +  Y2  - 

where  the  p.d.f.  of  is  given  by  d^  exp[-13y^/120],  -«•  <  y  <  °°  . 

Example  4.4.  Let  F  be  symmetric  Bernoulli  distribution  and  P  be  the  standard 
normal  probability  measure.  The  joint  distribution  in  this  example  is  given  by 

(4.30)  dC^U)  =  z^1(2ir)*n/2  [cosh(sn/n)  ]n  expC-yEx2]  Hdx . . 

In  this  example 

(4.31)  =  £°osh(z)}n> 

and 


(4.32) 


4/(z)  =  log  Ccosh(z)],  for  |  Re  a  1  ( z )  |  <®. 


The  large  deviation  rate  y(u)  is  finite  for  |u|  <  1  and  is  given  by 


(4.33)  y(u)  *  us2  -  log[cosh(s2) ], 

2 

where  tanh(s2)  =  u.  Since  the  c.g.f.  of  P  is  h(u)  =  u  /2 ,  we  get  for  u>0, 

(4.34)  V’(u)  =  s2  -  h'(<l>'(s2)) 

=  [s2  -  tanh(s2)]  >  0 

since  s2  =  tanh_1(u)  >  0.  It  is  easy  to  check  that  the  first  three  derivatives  of 
V  at  the  origin  are  zero  and  V^(0)  =  2.  Thus  zero  is  the  unique  global  minimum 
of  order  2.  This  completes  the  verification  of  Part  (I).  It  is  easy  to  verify 
that  the  lattice  random  variables  T^'s  satisfy  all  the  conditions  of  Theorem  3.16. 

2 

Part  (III).  Condition  (3.3)  is  trivially  satisfied  since  Cosh(x)  exp(-x  /2)  is 
an  integrable  function.  Also  for  £  >  0, 

(4.35)  l  exp[-£(V(u))]  5  ][  exp[£h(u)] 

u  e  R  ue  R 

n  n 

*  j  l  ^  exp[£u2/2] 

ue  R 
n 

1  ? 

*  5-  /  expC-tu/2]  -  0(n). 

-1 

Therefore  P  belongs  to  the  class  L*.  Thus  by  the  conclusion  of  Theorem  3.18  we 
get 

3/4  d  * 

(4.36)  Sn/n  ■”  Y2  » 

4 

where  V2  is  distributed  as  d2  exp(-y  /12),  .<*,  <  y  <  »  . 
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